A wide-field photometric survey for extratidal tails around five 
metal-poor globular clusters in the Galactic halo ^ 

Sang-Hyun Chun^, Jae-Woo Kim^, Sangmo T. Sohn^'^, Jang-Hyun Park^, Wonyong Han^, 
Ho-Il Kim^, Young- Wook Lee^'^, Myung Gyoon Lee^, Sang-Gak Lee^, and Young- Jong 

Sohn2 

Received ; accepted 

To be submitted to the Astronomical Journal 



^Department of Astronomy, Yonsei University, Seoul 120-749, Korea, sohnyj@yonsei.ac.kr 
^Institute for Computational Cosmology, Department of Physics, Durham University, 

South Road, Durham DHl 3LE, UK 

^Center for Space Astrophysics, Yonsei University, Seoul 120-749, Korea 
^Korea Astronomy and Space Science Institute, Daejeon 305-348, Korea 
^Astronomy Program Department of Physics and Astronomy, Seoul National University, 

Seoul 151-742, Korea 

''Current address: Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, 
MD 21218 



- 2 - 



ABSTRACT 



Wide-field deep gr'i images obtained with the Megacam of the Canada- 
Prance-Hawaii Telescope (CFHT) are used to investigate the spatial configura- 
tion of stars around five metal-poor globular clusters M15, M30, M53, NGC 5053, 
and NGC 5466, in a field-of-view ~ 3°. Applying a mask filtering algorithm to 
the color-magnitude diagrams of the observed stars, we sorted cluster's member 
star candidates that are used to examine the characteristics of the spatial stellar 
distribution surrounding the target clusters. The smoothed surface density maps 
and the overlaid isodensity contours indicate that all of the five metal-poor glob- 
ular clusters exhibit strong evidence of extratidal overdensity features over their 
tidal radii, in the form of extended tidal tails around the clusters. The orien- 
tations of the observed extratidal features show signatures of tidal tails tracing 
the clusters' orbits, inferred from their proper motions, and effects of dynamical 
interactions with the Galaxy. Our findings include detections of a tidal bridge- 
like feature and an envelope structure around the pair of globular clusters M53 
and NCC 5053. The observed radial surface density profiles of target clusters 
have a deviation from theoretical King models, for which the profiles show a 
break at 0.5 ~ 0.7rt, extending the overdensity features out to 1.5 ~ 2rt. Both 
radial surface density profiles for different angular sections and azimuthal num- 
ber density profiles confirm the overdensity features of tidal tails around the five 
metal-poor globular clusters. Our results add further observational evidence that 
the observed metal-poor halo globular clusters originate from an accreted satel- 
lite system, indicative of the merging scenario of the formation of the Galactic 
halo. 

Subject headings: Galaxy: halo — Galaxy: structure — globular clusters: general — 
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globular clusters: individual(M15, M30, M53, NGC 5053, NGC 5466) 
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INTRODUCTION 



In the current cold dark matter cosmology, the favorably accepted galaxy formation 
paradigm is that galaxies like the Milky Way were assembled from a series of accretion 



events involving lower mass subunits, in a fashion similar to ga . 



(IWhite fc Rees 



Ghigna et al 



1978 



1998 



Blumentah. 



Tormen et al. 



et a. 



1998 



1984 



Davis et al 



Ghigna et al 



1985 



axy cluster formatio n 



Tormen 



1997 



1998 



2000|). If these favored hierarchical 



models are correct, then the stellar halo of the Galaxy should form through tidal disruption 



and accretion of numerous low- mass fragments such as dwarf sate 



1998 




2005 





KlvDin et al. 


1999 


Abadi et al. 


2oodi 



Moore et al 



1999 



2006 



Bullock et al. 



lite galaxies (e.g., 



2001 



Johnston 



Bullock fc Johnston 



These merging events may leave traces of large-scale stellar 



substructures in the outer Galactic halo, and remnant stellar systems of the parent 
fragments with a long dynamical disruption time scale. 



Since the classical ISearle fc ZinnI (jl978l ) scenario of the formation of the Galaxy was 
first suggested, there has been more observational evidence to support the accretion scenario 
that the halo population of the Milky Way and globular clusters originated through the 
gradual merging and disruption of fragments. This evidence includes (1) a large-scale stellar 
substructure through the Galactic halo, (2) a feature of multiple populations in a globular 
cluster that is indicative of the remnant nuclei of a merged dwarf spheroidal galaxy, (3) an 
extended tidal feature in a spatial stellar density distribution, such as tidal tails and halos 
around a globular cluster, which might trace the cluster's orbit, and (4) a spatially aligned 
configuration of a group of several globular clusters. Consequently, globular clusters in the 



^Based on observations carried out at the Canada-France-Hawaii Telescope, operated by 
the National Research Council of Canada, the Centre National de la Recherche Scientifique 
de France, and the University of Hawaii. This is part of the Searching for the Galactic Halo 
project using the CFHT, organized by the Korea Astronomy and Space Science Institute. 
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Galactic halo and their surroundings, as the local manifestations of the hierarchical Galaxy 
formation, are the best places to look for the substructure of stellar distribution in the 
Galactic outer halo. 

For the first n ote of large-scale stellar substructures through the Galactic halo. 



Lynden-Belll ( 1l982l ) indicated that several dwarf spheroidal galaxies and outer halo clusters 
lie near two great circles passing close to the Galactic poles, i.e., the Magellanic stream 
and the Fornax-Leo-Sculptor stream. The spatially associated distribution of dwarf 
galaxies is interpreted as a result of the accretion of parent satellite galaxies which are 
considered the external fragments. Furthermore, recent large-sc ale surveys have uncover ed 



several complex substructures of the stellar Galactic halo (e.g., 



Ivezic et al 



2000 



Newberg et al 



200^ 



2009 



Martin et a. 



Yannv et al. 



all 1200. 
J 20061: 



2000, 



2003 



Rocha-Pinto et al 



2004 



Martinez-Delgado et al.l 



Juric et al. 



2005 



Grillmair &: Dionatos 



Belokurov et al. 



2006a 



20081 ) . The most striking example in this categor y is the discovery of the trails 



of ste 



2001a 



l ar debris associated with the Sagittarius dw arf galaxy (llbata et al. 



Majewski et al. 



2003 



Newberg et al. 



20031). 



1994 . 



1995 . 



19971 . 



Stellar population studies of globular clusters have suggested that the brightest globular 
clusters in the G alaxy are actually the remnant nuclei of dwarf spheroidal galaxies. Indeed, 



Lee et al 



(119991 ) detected multiple stellar populations of different age and metallicity 
in uj Centauri, the most massive globular cluster in the Galaxy. This is reminiscent 
of the Sagittarius dwarf galaxy, which includes M54 as its nucleus, the second most 



massive globular cluster in the Galaxy. Extending previous suggestions. 



van den Bergh 



( I2OOOI ) suggested that the majority of metal-poor young globular clusters in the outer 
halo originated from dwarf spheroida l galax ies, and it would be worthwhile to find their 



parent galaxies. Recently, iLee et al. 



(I2OO7I ) found that globular clusters with an extended 



horizontal-branch (EHB) would be kinematically decoupled from the normal clusters, and. 
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that these are the rehcs of the first building blocks of the Galaxy formation. 



If present Galactic globular clusters formed within larger stellar systems, they are 
likely surrounded by extratidal halos and tails made up of stars that were tidally stripped 
from their parent systems. Also, they would have lost the majority fraction of the initial 



ping 


Ostriker et al. 


1972 


Affuilar et al. 


1988; 


Kundic & Ostriker 


1995; 


Gnedin & Ostriker 


1997; 


Murali & Weinberg 


1997 




Gnedin et al. 


1999 




Baumeardt & Makino 


2003 


). This 



information suggests that surroundings around globular clusters provide an excellent 
example of such a structure, an d as a result we might expect to figure out them by 



wide- field deep photometry (e.g., 



Martinez-Delgado et al. 



Grillmair et al 



20041 ) . In t he first attempt to find 



3) 



( 119951 ) examined the stellar 



such a tidal feature around globular clusters, 
density distribution around 12 globular clusters and found that many globular clusters 
have surface density profiles significantly exceeding the prediction of formal King models 
and also extending outside the tidal radius by extratidal stars. Until now, signatures 



of the extratidal overdensit y feature have been found in over 30 globular c 



in the Galaxy and in M31 ( 



2000 



Testa et al. 



200G 



Grillmair et al. 



Siegel et al 



2001 



1996 



Lehmann fc Scholz 



Lee et al 



2003 



Sohn et al. 



.1997 
2OO3I : 



usters both 



Leon et al. 



Lauchner et al. 



20061 ). Recent star count surveys have revealed tidal structures stretching many degrees 



b e yond the cluster tidal rad i us, for which the most spectacu l ar examples include Paloma r 



5 (lOdenkirchen et al 



2001 



2003 



2009 



:lockosi et al 



2002 



Grillmair fc Dionatos 



20061) 



and NGC 5466 dBelokurov et al. 



2006b 



Grillmair fc Johnson 



20061 ) with tails extending 



as much as 22° and 45°, respectively. As a consequence of these previous studies, it has 
been apparent that tidal tails of globular clusters are likely aligned with the cluster's 
orbit because the debris produced by a dis solving globular cluster w ill have similar orbital 



parameters as the cluster itself. Recently, 



Odenkirchen et al. 



( 120091 ) presented kinematics 



of debris around Palomar 5 to suggest that the cluster's orbit in the sky is not exactly 
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aligned with its tidal tails. Here, we also note that iKoch et al.l (120041 ) found that the 
luminosity function (LF) in the tails of Palomar 5 agrees well with the LF in the outer 
region of the cluster but differs from the LF in the cluster's core. On the theoretical side 
many nu merical simulations 



clusters ( 



Dehnen et al. 



Kupper et al. 



2004 



lave reproducec 


a structure of 


;idal tails around globular 


di Matteo et al. 


2005; 


Lee et al. 


2006; 


Fellhauer et al. 


2007 



20081 ). Indeed, the cluster-galaxy interaction can cause stars to escape the 



cluster's gravitational potential, and then the unbounded stars remain in the vicinity of 



the cluster for several orbi 
behind the cluster's orbit 


;al periods to become tidal debris-forming tails in front of anc 


Combes et al. 


1999; 


Yim & Lee 


2002; 


CaDuzzo Dolcetta et al. 


2005; 


Montouri et al. 


2007 


)• 



The spatial alignment of the globular cluster system might also play an important role 
in untangling the possibility that the Galactic halo may have formed from the accretion of 
dwarf galaxies. For a statistically significant instance, the Sagittarius dwarf spheroidal has 
its own globular cluster s ystem associated with s everal halo clusters beyond a galactocentr ic 



distance Rqc > 10 kpc (iBellazzini et al. 



2003 



Tautvaisiene et al. 



2004 



Sb or done et al. 



20051 ). This attests to the fact that even m inor mergers 



and structure of the Galaxy. Furthermore, 



play a key role in the formation 



Yoon fc Led (120021 ) found that a group of the 



seven lowest metallicity ([Fe/H] < -2.0) clusters (i.e., NGC 5053, NGC 5466, M15, M30, 
M53, M68, and M92, classified as Oosterhoff group Il-b) displays a planer alignment in the 
outer halo perpendicular to the line joining the present position of the Sun and the Galactic 
center. They also suggested that these metal-poor clusters were originally born in the Large 
Magellanic Cloud (LMC) and have recently been captured by the Galaxy through the 
Magellanic plane. Indeed, the configuration of the Milky Way with the stream of clusters 
and satellites in the Magellanic plane closely resembles the stellar stream in the halo of 



M31, which lies along the sateUites M32 and NG C 205 (llbata 



low- metallicity clusters in the Magellanic plane (lYoon fc Lee 



et al 



2001bl ). Of the seven 



2002), NGC 5466 has been 
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the su bject of extensive investigations from wide-fie 



data (IBelokurov et al. 



2006b 



Grillmair fc Johnson 



d ph otometric analyses of the SDSS 



20061 ) to show a feature of extended 



extratidal tails outside the tidal radius along the cluster's orbit. The other clusters in the 
Magellanic plane have not yet shown any strong observational evidences of extratidal tails 



19951), M68 flGrillmair et al. 



19951), and M92 flTesta et al. 



2000 



1 

.. M15 (iGrillmair et al. 


Lee et al. 


2003 


)• 



In this paper, we investigate the two-dimensional distribution of stars around five 
clusters (M15, M30, M53, NGC 5053, and NG C 5466) of the sev en extreme metal-poor 



2OO2I ) and we report the spatial 



globular clusters of the Oosterhoof group Il-b (lYoon fc Led I: 
structure of the stellar distribution in the vicinity of the clusters. Wide-field g'r'i' imaging 
photometric data for the five globular clusters have been secured from the observations of 
the Megacam attached to the Canada- France-Hawaii Telescope (CFHT). Table [1] contains 
basic data for the five target globular clusters. Note that the other two low-metallicity 
clusters M68 and M92 have not been observed at the present observing runs. Section 2 
presents the observations, data reduction procedure, and photometric measurements of 
the resolved stars. The statistical mask filtering method for the photometric selection of 
clusters' member stars is presented in Section 3. In Section 4, we investigate extratidal 
features around the target clusters by using two-dimensional stellar density maps, radial 
surface density profiles, and azimuthal number density profiles. The results are discussed 
and summarized in Section 5. 



2. OBSERVATIONS AND DATA REDUCTION 



The wide-field photometric imaging data for M15, M30, M53, NGC 5053, and NGC 
5466 were taken in Queued Service Observing (QSO) of the 3.6 m CFHT. The detector in 
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the Megacam is a mosaic of thirty-six 2048 x 4612 pixels^ CCDs, that together cover an 
area of ~ 1° x 1° wide-field-of-view with a plate scale of 0.187" pixel"^. We covered nine 
and eight Megacam observing fields for the surroundings of two clusters, M15 and NGC 
5466, and of one cluster, M30, respectively. For each cluster, one field was centered to 
imbed the cluster itself, while the other fields outside the object were chosen so as to cover 
a total field of view of ~ 3° x 3°. Since the tidal radii of the clusters are in the range of 14' 
~ 22 , a set of these observations provided sufficiently large spatial coverage for our purpose 
of studying the stellar density distribution and tidal features in the vicinity of the globular 
clusters. The other two clusters, M53 and NGC 5053, compose a close neighbor pairing 
system with an angular separation of ~ 1° on the sky, and the two clusters have nearly the 
same distance modulus. Note that the close proximity of the two clusters is very unusual in 
the outer halo of the Galaxy. In order to search for direct observational evidence of a tidal 
interaction between the two clusters and tidal features around each cluster, we observed 
nine Megacam fields in a pattern as shown in Figure [H covering both of the clusters. 

The imaging data for the target clusters were taken in three SDSS filters, g ^ r', and 
i', respectively. Single exposures were taken for each filter without a dither pattern. The 
exposure times to be used were determined so as to achieve a limiting magnitude of ~ 3 
magnitudes below the main sequence turnoff of each cluster, i.e., g ~ 23 mag. Table [2] 
presents the journal of the observations. Stars have 0.6 — 1.0 arcsec FWHM in the final 
processed images, depending on the filter. 



The removal of instrumental s ignatures from the raw 



processing pipeline at the CFHT (IMagnier fc Cuillandre 



image s was done with the ELIXIR 



20041 ). This processing includes 



bias subtraction, fiat fielding, and the subtraction of a fringe frame. Also, the pre-analysis 
of the ELIXIR produced 36 CCD images of a given mosaic with the same zeropoint and 
magnitude scale. A typical value of the residuals across a mosaic is ~ 0.0086 mag. 
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The photometric measurements were made with the 



poin t 



1987 



-spread-function 



Stetson fc Harris 



(PSF) 



19881 ) 



fitting programs DAOPHOT II and ALLSTAR fiStetsonI 
Analysis was done on a single-chip basis for all exposures of the 36 CCD chips on the 
mosaic. A PSF varying quadratically through positions was constructed using ~ 100 — 200 
bright and isolated stars for a heavily crowded image and ~ 25 — 50 stars for a less 
crowded image. The quality of each PSF has been improved by removing neighboring stars 
around selected PSF stars and reconstructing the PSF iteratively. Instrumental magnitudes 
for each star based on the determined PSF were then measured by the ALLSTAR 
process. Ape rture correction was applied by using the growth-curve analysis package 



DAOGROW ( Stetson fc Harris 



1988 



Stetson 



19901 ) to compensate for the missing light 



lying outside of the PSF tail for ev ery star in ea ch frame. After the aperture correction 



DAOMATCH and DAOMASTER fIStetson 



19921 ) were used to match stars on each CCD 



chip, and coordinate offsets for each field were estimated from the positional data in each 
frame. We then converted all the positional data of each chip on a mosaic into a single 
coordinate system. Finally, the positional data for stars in the observed Megacam fields 
were transformed i nto the equatorial coordinate system, by using the USNO-B catalog 



( iMonet et al. 



20031 ) for bright stars in the fields. 



The measured instrumental magnitudes were photometrically calibrated using the 
zeropoints that were computed from standard star observations during each Megacam run. 
The calibrated magnitudes were fi nally converted into a standard SDSS photometric system. 



as defined by 



Smith et al. 



( 120021 ). Individual stars with poorly determined photometric 
measurements, as reported by ALLSTAR, were removed to limit the number of spurious 
detections for each cluster. Objects with x ^i^nd SHARP values that markedly deviated 
from the values for the majority of objects were rejected, and this removed objects near 
the faint limit of the data and extended objects such as faint galaxies. In order to correct 



the interstellar extinction, we derived the individual extinction value for each detected star 
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from the maps of ISchlegel et al.l (Il998l ). and these values have been subtracted from the 
observed magnitude. Finally a total of 239 600, 100 069, and 69 901 stars were detected for 
each ~ 3° X 3" field of clusters M15, M30, and NGC 5466, respectively. For nine Megacam 



fields of M53 and NGC 5053, we detected 79 075 stars in all three g r i images. 



PHOTOMETRIC SELECTION OF THE CLUSTER STARS 



A feature of spatial stellar distribution around globular clusters represents the sum 
of the unbound stars released by the cluster and the field stars. Therefore, to trace the 
overdensity features around the clusters in detail, one needs to minimize the number of 
contaminating field stars in a given sample and to enhance the contrast b etween the clust er 



and t 



Grillmair et al. 



le field. In the first statistical attempt to discriminate field stars, 
( ll995l ) introduced a polygonal mask method in a color-magnitude (C-M) space in which 
only field stars with colors and magnitudes that resemble those of cluster stars are counted 
to optimize the ratio of the cluster stars to field stars. However, the defined mask in a C-M 
space does not provide the optimal method to map the spatial stellar distribution around 
sparsely populated globular clusters because many bins with a small number of stars in a 
C-M space can be statistically exclu ded from the mask. To reduce the disadvantage of poor 



statistics, 



Odenkirchen et al. 



(1200 ll ) introduced a new method to compress the photometric 



Odenkirchen et al. 



fl200lh to 



data in a C-M space by constructing new color indices from which the weighted sum of field 
stars plus cluster stars yields the best approximation of the observed total C-M distribution. 
In the present study, we have applied the same method of 
define a C-M mask with new color indices C\ and C2 and g r % magnitudes for detected stars 
in the vicinity of each cluster. 

First, we defined two orthogonal color indicies C\ and from a linear distribution of 
stars in a i^g — r\ r' — i') color-color diagram. Indeed, the choice of the indicies is such 
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that the main axis of the almost one-dimensional locus of cluster stars in the color-color 
plane lies a long the C\ axis, while the axis is perpendicular to it (e.g. 



Odenkirchen et al. 



2001 



20031 ) . Equation ([T]) presents general forms of the rotational transformation of the 



((?' — r' , r' — i') plane to the coordinate system of (ci, C2). Here, the systematic variation of 
color with magnitude is mostly contained in c^ , while variations in ar e mostly due to the 



observational photometric errors, as noted in I Odenkirchen et al. 



(120011 ). Transformation 



coefficients a and b for each cluster were determined by using the least-squared fit to the 
data, which are listed in Table 3. Assuming that two clusters M53 and NGC 5053 have 
very similar color-magnitude diagrams(CMDs) w ith almost identical metallicity and age 



[e.g. 



Heasley fc Christian 



1991 



Rey et al. 



19981 ). and nearly the same distance modulus. 



we estimated a set of Ci, C2 and a, b for all stars detected in both of the clusters. 



Ci = CL{g — T ) + b{r — i ), 
C2 = —b{g — r ) + a(r — i') 



(1) 



In (c2,/) space of stars located in the central region of the clusters, we preselected the 
sample in C2 by discarding all stars with \c2\ > 2(t c^ii)., where (Jc2{'^') is the photometric 
errors in C2 for stars with magnitude i' . The radii of the cluster's central regions were 
empirically defined in the unit of half-mass radius (r/^) or core radius (r^) of the cluster in 
order to find appropriate color-magnitude sequences of stars in each cluster. The heavy 
lines in the leftmost panels in Figure 2 give 2(Tc2(«') limits for stars in each cluster. For a 
system of M53 and NGC 5053, we selected a central region of the clusters as the sum of 
3rh<r < Svh for M53 and < r < for NGC 5053. 

For s tars preselected in th e {c2,i') plane, we then applied the C-M mask filtering 



method (iGrillmair et al 



19951 ) in the (ci,z') plane, from which we eventually extracted a 



sample of stars to trace their spatial configuration around the clusters. We first examined 
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the distribution of the cluster population in the (ci, i') space by means of preselected stars 
in the cluster central area and the background field areas. The cluster central area was 
assigned so as to have the largest number ratio of stars enclosed 2o"ci(z') and those outside 
2ac^{i) in the (ci,/) plane for all preselected stars within the radius of the circle. Here, 
crci(i') is an estimated standard deviation of i magnitudes of stars at ci. The second from 
the left panels in Figure [2] show (ci,i') CMDs of stars within the radii of the determined 
central areas of the clusters, which are denoted by circles centered on each cluster in 
Figure [31 We also selected four background field areas in circles at sufficiently distant parts 
from the cluster. For three isolated target clusters, M15, M30, and NGC 5466, the sum of 
the background field areas is nine times those of the determined cluster's shown 
in Figure [3l Note that, as shown in the bottom-right panel of Figure [3l the background 
areas for a field of M53 and NGC 5053 were empirically selected to have a total area five 
times larger than the sum of central areas of the two clusters. The third and forth panels of 
Figure [21 from left to right, show the (ci,z') CMDs of the stars in the selective background 
regions and in the total survey region for each cluster, respectively. 

Next, we compared the (ci,i') CMD of stars in the selective cluster central region 
with that of stars in the assigned background regions, from which we select areas in the 
(ci,z') plane with high significance about the expected true number of cluster stars. This 
procedure was done by subdividing the C-M plane into small subgrid elements with a 
0.025-mag width in ci and a 0.12-mag height in i' . Then, we measured the number of stars 
in each subgrid element of cluster area nci{ci,i') and those of background fields n/(ci,i'), 
for which we assumed that the distribution of background field stars on the C-M plane does 
not vary across the observed field. The local signal-to-noise ratio s(ci,/) in each subgrid 
element was then calculated by Equation ([2]), where g is simply the ratio of the area of the 
cluster region to that of the background field. 
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s{Ci,t j = (2) 
^/nci{ci,i') +5(%/(ci,z') 

From the array s, we obtained an optimal mask envelope in the (ci,i') plane by setting 
a threshold < Smax (maximum of s) and identifying the subgrid elements of s > sum- 
In order to do this, the elements of s{ci,i') were sorted into a series of gradually decreasing 
thresholds over the one-dimensional index k. From the determined subgrid element with 
the highest signal-to-noise ratio (smax), a cumulative number of stars in the cluster area 
Ndia-k) = X]f=i^c/(0 '^as then counted through a progressively larger area = fca;, where 
ai is the area of a single element of the C-M plane. The cumulative star counts for the 
background field Nf{ak) = Yl^=i "^fO) were calculated in the same manner as for the cluster 
area. The cumulative signal-to-noise ratio S{ak) was then determined by Equation ([3]): 



Sia ) = ^^'(^fe) ~ aNfjak) 

S{ak) reaches a maximum value for a particular subarea of the C-M plane, and the 
corresponding value of s{ci,i') is set to a threshold, i.e., sum- Now, the filtering mask area 
is chosen by selecting subgrid elements with larger s(ci,i') values than the determined sum- 
The heavy lines in (ci,i') CMDs in Figure [2] represent the boundaries of the selected subgrid 
elements. The faint limits of the filtering mask areas were assig ned to be i = 22.0 for NGC 
5466 and i' = 22.5 for the other clusters to avoid biases caused by the poor completeness 
(~ 70%) of the photometry and by the variable seeing conditions over the surveyed sky 
areas during the observations. Finally, the entire sample of stars in the determined filtering 
mask area was considered in the subsequent analyses to examine the spatial configuration 
around the target clusters. 
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4. SPATIAL CONFIGURATION OF STARS 

In this section, we present and discuss the analytic results related to the spatial 
configuration of stars in the vicinity of target clusters individually for each cluster. Details 
of the extended tidal features with significant overdensities of stars over the limiting tidal 
radii of the clusters are examined by (1) the spatial surface density map of stars, (2) the 
radial profiles of the surface density, and (3) the azimuthal number density profiles as a 
function of position angle at various radial distances from the cluster center. 

To construct a spatial surface density map, we first made a two-dimensional star 
count map by binning into small grids on the spatial stellar distributions of the C-M 
mask-selected stars for each cluster. The background density level was simply assigned 
by the measured mean surface density value for the region outside the tidal radius of 
each cluster since a variation with relatively low background density does not seriously 
affect to the two-dimensional configuration of stellar distribution around a cluster. The 
background-subtracted star count map was transformed into a smoothed surface density 
map by means of a Gaussian smoothing algorithm that enhanced the low spatial frequency 
of the background variation and removed the high frequency of the spatial stellar density 
variation. Then, we overlaid isodensity contour levels with a standard deviation unit (a) of 
the background level on the smoothed maps with various kernel values. 

As an analytic indicator to confirm the overdensity feature around the tidal radius 
of a cluster, the radial surface density profile of a cluster was derived by measuring the 
number of stars in concentric annuli with a l.O' width. The effective radius of each annulus 
is assigned by the equation r^, = A/0^5(rf~Tr|^, where and rj+i are the inner and 
outer radii of an annulus, respectively. Radial surface density profiles for a different 
direction are also derived, for which we divided an annulus with 3.0' widths into eight 
sections (SI— S8) with an angle of 45°, as shown in Figure HI Here, we caution that the 
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radial surface densities at the central region of cluster fields are usually underestimated by 
crowding effects. In order to examine the crowding effects, we added appropriate number 
of artificial stars with i magnitudes in the range of the C-M mask-selected area (Figure [2]) 
to the i images of each cluster, and we calculated the mean radial recovery rates of the 
input artificial stars as completenesses in photometric measurements out to ~ 25' from the 
clust er's c e nter. Considering the determined photometric completenesses, we empirically fit 



the 



KingI (119661 ) models to the observed radial surface density profiles of target clusters 



over a radial range t hat keeps up recovery rates in stable with high com pletenesses. It 



is worth noting that 



McLaughlin fc van der Marell (120051 ) have used the IWilsonI (119751 ) 



models to describe properly the outer surface density structure of the majority of globular 
clusters in the Galaxy and the Ma gellan i c Clo uds. It is true that the Wilson models are 



spatially extended more than the 



King] (119661 ) models, but the Wilson models have still 



finite truncated tidal radii. This leads that the Wilson models are unlikely to describe 
extratidal overdensity features of globular clusters, which are usually characterized by a 
break and radially untruncated power laws of the surface density profiles. 

Clusters with obvious extratidal extensions showed a break in the observed radial 
surface density profiles near the tidal radius ( 



Testa et al. 



2000 



Rockosi et al 



from the form predicted by the 



2002 



^ee et al 



Grillmair et al 



20081: 



1995 



Olszewski et al. 



;eon et al. 



2000 



20091 ). which departs 



KingI (119661 ) models. Numerical simulations also confirmed 

le features 



such a break with power laws in the ra dial surface density profiles due to ; 
of tidal debris around the clus ters (e.g. 



Indeed, 



Johnston et al. 



Combes et al. 



1999 



Johnston et al. 



1999 



2002 ). 



(119991 ) predicted that for a constant mass-loss rate, the radial 



surface densit i es of stripped stars in the tidal tails are proportional to r ^. However, 



Rockosi et al.l (l2002l ) found that the annular-averaged density of stars along the well-defined 



extratidal tail s around Pal 5 i s fit to a power law with a significantly steeper slope than that 



predicted by 



Johnston et al. 



(119991 ). For the target clusters in this paper, we fit a power 
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law of r"^ to the outer part of radial pr ofiles measured in fu ll conc entric annuli, and com pare 
the measured slope (7) with those of 



Rockosi et al 



(120021) and 



Johnston et al. 



fll999f ). To 



investigate directions of the overdensity feature, we also measure slopes of power laws for 
the radial surface density profiles to eight angular sections with different directions. 

The azimuthal number density of stars with respect to a position angle was also derived 
for stars in an assigned annulus. The position angle was measured clockwise from the east 
principal axis in a 10° unit. Widths of each annulus were assigned to be lO' and 15' for the 
inner and outer regions, respectively, in order to contain a statistically sufficient number 
of stars in each area. The radially cumulative azimuthal number density profile was also 
derived to trace the direction of the stellar overdensity feature around the clusters. 



4.1. M15 

Figure shows a raw star count map around M15 and the surface density maps 
smoothed with various Gaussian kernel values of 0.09°, 0.135°, and 0.195°, from the 
upper-left panel to the lower-right panel. The spatial bin size of the star count map was 
originally set to be 1.8' x 1.8'. Several individual blank chips, which did not contain useful 
observational data in Figure O were excluded from the subsequent analyses. Isodensity 
plots were overlaid on the maps with the contour levels of 2cr, 2.5cr, 3a, 4cr, 80", and 30cr above 
background density level, where a is the standard deviation of the mean background value. 
Contours in the star count map are the same as those in the smoothed surface density 
map with a Gaussian kernel value of 0.195°. The different arrows indicate directions of the 



/irv co s 6 = —0.30 ± 1.00 mas yr ^ and fis = —4.20 ± 1.00 mas yr ^ ( ICudworth fc Hanson 
1993 ). The circle centered on the cluster indicates its tidal radii 
IBM). 



(Cudworth & Hanson 


f n = 21.5' ( 


Harris 



- 18 - 



It is apparent in Figure [5] that the overdensity feature is obviously extended over the 
tidal radius of M15 at the levels larger than 1o above the background. Furthermore, the 
extratidal overdensity feature beyond the tidal radius seem to form extended tails around 
the cluster. Among them, long tails to the southwest and northeast directions around the 
cluster are aligned apparently with the directions of the Galactic center and anticenter. 
Marginal extensions are present in the outskirts to the south and north directions, which 
could correspond to the cluster's orbit. The overdensity feature of the tails becomes more 
prominent in the lower resolution surface density maps with a larger Gaussian kernal value. 
It is noted here that we did not find distinctive evidence for a southeast tidal extension as 
observed by 



Grillmair et al. 



(119951 ). The difference could be accounted for by the lower 



spatial resolution and the brighter photometric limiting magnitude used b y them. Indeed, 



most of the cluster stars in this study are faint main-sequence stars, while 



Grillmair et al. 



( 119951 ) used the subgi ant branch sta r s and main-sequence stars near the turn-off point of 



the CMD. Note that 



Combes et al 



(119991 ) certified in their simulation that low mass stars 



were usually stripped from a globular cluster and formed tidal tails or streams. 

The radial surface densities for M15 measured in each concentric annulus are 
plotted in the upper panel of Figure [6] along with a theoretical King model of c = 2.5, 
which are arbitrarily normalized to our measurements. We also show two plots of the 
observed-minus-predicted residuals, log{0/C), to illustrate the departures of the observed 
radial surface densities from the prediction of the King model, and of the photometric 
recovery rates to consider radial crowding effects on the observed profile. It is apparent in 
the plots that incompleteness due to crowding effects reduces the measured surface densities 
in the innermost region, while the radial surface density profile in the region outer than 
log{r') ~ 0.7 has completeness values larger than ~ 95%. Consequently, the radial surface 
density profile departs from the King model with a break at a radius less than the tidal 
radius, i.e., ~ O.Tr^, and the excess extends to ~ 1.5rt. The radial surface density profile 
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in this overdensity region resembles a radial power law with a slope of 7 = — 1.5 9 ± 0.21 



which is comparable to the value —1.58 ± 0.07 proposed by 



Rockosi et al 



fl2002h from the 



annular-averaged densities of stars along the extratidal tails of Pal 5. However, the slope 
is significantly ste eper than the case o: 



over a long time (jJohnston et al. 



7 = — 1, predicted for a constant mass-loss rate 
I999I ). Radial surface density profiles to eight angular 
sections with a different direction, as illustrated in Figure HI are also represented in the 
lower panel of Figure [61 The overdensity feature in the region of O.Tr^ ^ r < l.Sr^ is 
commonly detected in all eight angular sections. The estimated mean surface densities (/i) 
in the angular sections 1, 4 and 8 are higher than those in the other sections. This is in 
good agreement with the extended overdensity features to the directions of the Galactic 
center and anticenter, as shown in the surface density maps of Figure [51 Comparing the 
power laws of the radial surface density profile in each angular section, the slopes (7) seem 
to be shallower in angular sections 2 and 6, which are likely aligned with the orbit of M15, 
i.e., the direction of the cluster's proper motion. 

Figure [7] shows the annular and cumulative azimuthal number density profiles around 
M15. As presented in Figure [H we measured the number density clockwise from the east 
principal axis with a 10° angular bin. Note that number densities at the angular bins 
without data, caused by the blank chips, are interpolated by the values of adjacent angular 
bins. The radial range 15' ^ r < 75' corresponds to O.Tr^ < r < 3.5rf, in which we examine 
the direction and the extension of the tail feature of M15. The radial bins are set to be 
10' for the inner three annuli and 15' for the outer two annuli. At the innermost annulus 
with 15' ^ r < 25', it is apparent that the number densities in the ranges of position angles 
~ 150° to ~ 220° and ~ 300° to ~ 60° are clearly higher than those in the other position 
angles. The high density features in the positions are likely to extend to the whole observed 
field on the sky, as shown in the cumulative azimuthal number density profile. This is in 
accordance with the tail features to the direction of the Galactic center and anticenter, as 
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shown in the surface density map of Figure [51 With a somewhat low confidence, there also 
seems to be relatively weak overdensity features at position angles ~ 100° and ~ 280° in 
the regions of the inner two annuli (i.e., 15' ^ r < 35'). These are reliably matched with 
the tails to the direction along the cluster's orbit, as inferred from the surface density maps 
and the radial density profiles of M15. 



4.2. M30 

In Figure [8] we show the spatial configuration of stars surrounding M30. The top-left 
panel shows the two-dimensional star count map that was made by measuring the number 
of cluster stars in each 2.5' x 2.5' bin on the C-M mask-selected star count map of M30. The 
other panels to the bottom-right show the smoothed surface density maps with Gaussian 
kernel widths of 0.084°, 0.147°, and 0.210°. Note that in Figure [8]the Megacam observations 
were not secured at the outer southeast field of M30. The isodensity contours with the 
background level and 0.5cr, 1.5a, 2.5cr, 3.5a and 8a levels above the background were overlaid 
on the surface density maps. Contour levels in the star count map correspond to those 
of a smoothed map with a Gaussian kernel width of 0.210°. The short arrow indicates 
the projected direction to the Galactic center on the sky, and the long arrow indicates 
the projected proper m otion with /x^ cos S = 1.42 ± 0.69 mas yr~^ and fis = —7.71 ± 0.65 



fll999h . The tidal radius of M30 (i.e., n = 18.34') 



mas yr ^, ineasure d by iDinescu et al. 

from iHarrisI (119961 ) was also plotted as a circle centered on the cluster. 



At contour levels larger than 0.5cr above the background, it is apparent in Figure [8] 
that a flocculent overdensity feature with tails is well developed around M30 and extends 
to about 2rt ^ 40' and beyond. While there is no strong evidence of extratidal features 
extending into specific directions in the surface density maps, the apparent east-west 
extension seems to be aligned with the direction of the Galactic center and anticenter. 
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There is no apparent extended extratidal tail toward the direction of the proper motion on 
the cluster's orbit in the surface density maps. Instead, distant stars unbound from the 
cluster are likely building up a feature of the northern extension to the opposite direction 
of the cluster's proper motion, with significant overdensities up to a contour level of 2.5o" 
in Figure [HI The northern extension feature is likely tracing the cluster's trailing tail 
through its orbital path. We note, however, given the position and distance of M30 from 
the Galactic plane, the northern ex tension of the star s could be a result of gravitational 



interaction with the Galactic plane (jPiotto et al. 



19971) 



A radial surface density profile for M30 was derived by measuring the mean number 
density of stars in concentric annuli and is presented in the upper panel of Figure [H For 
the inner profile, we overplotted a theoretical King profile of c = 2.5, which was arbitrarily 
normalized to the observed radial surface density profile of M30. The observed- minus- 
predicted residuals illustrate the deviations of the observed radial surface densities from 
the prediction of the King model. The innermost region shows crowding effects out to 
log{r') ~ 0.5, after which the photometric completenesses keep a nearly constant value of 
~ 98%. A break in the observed radial density profile appears at a radius less than the tidal 
radius of M30, i.e., ~ 0.6rf, and the overdensity feature extends to about 2rt of the cluster. 
The departure of the observed radial surface density profile in this region from the King 
profile is characterized by a power law with a slope of 7 = —1.41 ± 0.21. Similar to the case 
of M15, the slope is also significantly steeper than the inde x of predicted for extratidal 



structure with a constant mass- loss rate over a long time (1 Johnston et al.l 119991 ). Radial 
surface density profiles to eight angular sections are plotted in the lower panel of Figure [9l 
Note that incompleteness due to crowding reduces the measured surface densities at the 
innermost point with respect to those expected by the King model. It is apparent that the 
profiles have in common a break at ~ O.Gr^ and an extended overdensity feature to ~ 2rt 
with a power law. At the radial range, the estimated mean surface densities are relatively 
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higher in the angular sections 4 and 8 than in the other sections. This can be explained 
by the spatially extended extratidal features to the directions of the Galactic center and 
anticenter, as shown in the surface density map of Figure [Hi We also note that the slopes of 
the radial density profiles of the angular sections 2 and 6 are shallower in the radial range 
of 0.6rt < r < 2rj than those of the other angular sections. This feature may be related to 
a cluster's orbit inferred from the proper motion, although the spatial surface density maps 
do not show an extinctive extratidal feature of stars on the cluster's orbit. 

The left and right panels of Figure [TD] show the annular and radially cumulative 
azimuthal number density profiles, respectively, of M30 with an angular bin of 10°. The 
radial bins are set to be lO' for the inner three annuli and 15' for the outer two annuli. 
The tidal radius of M30 is placed in the innermost annulus. As shown in Figure [10] the 
azimuthal number density profiles have a flocculent feature with respect to the position 
angle, which is also expected in the spatial surface density maps of Figure [81 Nevertheless, 
there appear to be two main overdensity ranges of the position angle, i.e., ~ 120° to ~ 220° 
and ~ 300° to ~ 40°, in the azimuthal number density profile for the innermost annulus. 
Also, the overdensity feature in the ranges of the position angle extends to the second 
annulus, i.e. < 2r(, with a somewhat low confidence. This might be associated with the 
apparent east-west extension of stars to the direction of the Galactic center and anticenter, 
as shown in the spatial surface density maps in Figure [H] In the azimuthal number density 
profiles, we did not find evidence of the extended extratidal features induced by the cluster's 
orbital motion at the remote area in the vicinity of M30. 



4.3. M53 and NGC 5053 



Figure [TT] represents the spatial configuration of selected stars in the vicinity of a pair 
of the two clusters M53 and NGC 5053. As described in Sections 2 and 3, cluster stars for 
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both clusters are selected by using the same selection criteria in C-M spaces of (ci,z') and 
(c2,i'). The upper-left panel is a raw star count map with a spatial bin size of 3.0' x 3.0'. 
The surface density maps smoothed with various Gaussian kernel values of 0.125°, 0.175°, 
and 0.240° are also plotted in Figure [11] to bottom-right. The long arrow indicates the 
direction of the proper motion of M53, and the two short arrows represent the direction 
to the Galactic center. The proper motions of M53, i.e., UnC o sS = 0.50 ± 1.00 mas yr 



/i5 = —0.10 ± 1.00 mas yr ^ are adopted from 



Dinescu et al 



(119991 ). The proper motion of 



NGC 5053 has not yet been reported. C ircles are tid al radii for each cluster, i.e., = 21.75 



for M53 and n = 13.67' for NGC 5053 (IHarris 



19961). 



It is remarkable in Figure [TT] that the spatial distribution of stars in the vicinity of 
M53 and NGC 5053 show a complex configuration with various extratidal features. The 
contour levels are background level, and 2a, 4ct, 8a, and 20a above background level, where 
a is the standard deviation of the estimated background level. The main new findings that 
were observed in the surface density maps are as follows: (1) the extratidal overdensity 
features extended over the tidal radii of both clusters at the contour levels above the 
background, (2) features of tails on the cluster's orbit in the east-west direction for M53 
and to the direction of the Galactic center for NGC 5053, (3) complex substructures of 
extratidal feature of clumps and ripples around the clusters, (4) a tidal bridge-like structure 
connecting spatially the close neighbor pairing of M53 and NGC 5053, and finally (5) an 
envelope composed by stars surrounding the cluster system. Of particular interest are 
the tidal features that are possibly enhanced by the dynamical interaction between two 
clusters, which include features such as clumps and ripples surrounding the clusters and 
most strikingly, a tidal bridge-like feature and an envelope structure around the multiple 
globular clu ster system. Note th at there are no known binary globular clusters in the Milky 



Way, while 



Minniti et al. 



(j200J) reported the discovery of a binary cluster in the peculiar 



giant elliptical galaxy NGC 5128. We therefore point out that these observed extratidal 



-24- 



structures of the pairing system of M53 and NGC 5053 could indeed play a unique key 
role in the effect of near passage on the dynamical tidal interaction between two globular 
clusters in the Galaxy. 

The upper panel of Figure [12] shows the radial surface density profile of M53, which was 
derived by measuring mean stellar number densities in each concentric annulus. Stars within 
the tidal radius of the companion cluster NGC 5053 were not included in the radial surface 
density profile of M53 to avoid overlapping of stars in NGC 5053. In order to compare the 
observed radial density profile, a theoretical King model was computed with a concentration 
parameter c = 1.78. As apparent in Figure [T21 the observed profile departs from the King 
model with a break at a radius less than the tidal radius, i.e., 0.7rj. Note that the surface 
densities in the innermost region is affected by crowding effects out to log(r') ~ 0.8, after 
which the photometric completenesses show a nearly constant value of ~ 98%, as shown in 
the plot of the radial photometric recovery rates. The overdensity feature, which extends to 
l.Qrt, is likely represented by a power law with a slope of 7 = —1.58 ± 0.19, which is almost 



identical with the value —1.58 ± 0.07 



extratidal tails of Pal 5 ( Rockosi et al. 



' or the annular-averaged densities of stars in the 
2OO2I ). but significantly steep er than that predicted 



I999I ). The residual plot 



from a constant orbit-averaged mass loss rate (iJohnston et al. 
within the tidal radius also exhibits a departure of the observed radial surface densities 
from the predic tion of the theoretica l King model. As a signature of the presence of tidally 



stripped stars. 



Beccari et. al. 



(I2OO8I ) also recognized a change in slope of the radial surface 



density profile for M53. 

The lower panel of Figure [12] shows radial surface density profiles of M53 in different 
directions. Of the total eight angular sections defined in Figure H] we could not derive the 
radial profiles of the angular sections 5 and 6 because observational data were not secured 
for the southern area of M53. As can be seen from Figure [12], the obtained radial profiles for 
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the six angular sections show behaviors of a break at ~ 0.7rj in common, becoming power 
laws over the tidal radius to l.Gr^. This trend resembles the extended overdensity features 
over the tidal radius appearing in the surface density maps of Figure [TTJ The estimated 
mean densities (/i) at the radial range show higher values in the angular sections 1, 4, and 
8. This feature corresponds to the apparent extratidal tails on the cluster's orbit in the 
east-west direction aound M53, as shown in Figure [TTl Also, slopes (7) of the power laws 
in the angular sections 3 and 4 seem to be slightly flatter than the mean slope of the radial 
surface density profile of M53. It can be noted that angular sections 3 and 4 are equivalent 
in the directions of the Galactic anticenter and the trailing tail on the cluster's orbit. 

In Figure [131 we show the annular and cumulative azimuthal number density profiles 
around M53 with a 10° angular bin. We have divided the total sample in the radial range 
of 15' ^ r ^ 75', which corresponds to O.Tr^ < r < 3.5ri, into five concentric annuli with 
radial bins with lO' and 15' widths for the inner three and outer two bins, respectively. 
The tidal radius of M53 is placed in the innermost annulus. The blank bins in Figure [13] 
correspond to the areas without photometric data and the area within the tidal radius of 
the neighbor cluster NGC 5053. The azimuthal density profiles for the two inner annuli 
indicate the overdensity feature in the range of position angles ~ 160° to ~ 240°, which 
is in good agreement with the direction of the trailing tail on the cluster's orbit. Also, an 
apparent overdensity feature in the range of position angles ~ 310° to ~ 20° likely traces 
the tidal bridge-like feature in the direction of the companion cluster NGC 5053, as shown 
in the surface density maps of Figure [HI 

We now turn to NGC 5053, w hich is a very diffuse cluster with a low concentration 



parameter c = 0.85 (IHarris 



1996! ) . Figure [TH shows the derived radial surface density 
profiles of NGC 5053. As in the case of the companion cluster M53, stars within the tidal 
radius of M53 were not used to derive the radial density profile of NGC 5053. The upper 
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panel shows the profile of the mean surface densities measured in each concentric annulus, 
and the lower panel consists of those of eight angular sections. The observed profiles were 
compared with a theoretical King model of c = 0.85. By crowding effects, surface densities 
seem to be reduced in the central region of the cluster inner than log{r') ~ 0.5, after which 
the photometric completenesses show a nearly constant value of ~ 96%. A break at ~ O.Tr^ 
is commonly appeared in the obtained radial density profiles. The overdensity feature 
extended to 2.5rf of NGC 5053. For the radial profiles measured in concentric annulii, the 
departure from the King model in this region can be represented by a power law with a 
slope 7 = —0.62 ± 0.15. It is remarkable that the radial surface density profiles of NGC 
5053 are characterized by a long extension of the extratidal overdensity feature and a much 



shallower slope o 
mass-loss rate of 



: the power law compa red with the case of the slope 7 = —1 in a constant 



Johnston et al. 



of the cluster with a small concentration of c = 0.84 (iHarri, 



( 19991 ). This is indeed consistent w ith the diffuse structure 



19961 ). In the overdensity 



region, the mean surface densities [fi) are higher in the angular sections 3 and 4, and the 
angular sections 7 and 8 than in the other angular sections. Angular sections 3 and 4 are 
coincident with the tidal bridge-like feature toward the neighbor cluster M53. Also, the high 
density in angular sections 7 and 8 might be related to the apparently extended extratidal 
features in the direction of the Galactic center, as shown in Figure [Til 

Figure [15] shows the annular and radially cumulative azimuthal number density profiles 
around NGC 5053, as a function of position angle. The radial range 10 < r < 80 
corresponds to 0.7r( < r < 5.9rt, in which we set six concentric annuli. The blank bins 
correspond to the areas without photometric data and the area within the tidal radius 
of the neighbor cluster M53. It is apparent that the azimuthal number density profile 
shows a flocculent feature at the innermost annulus in which the tidal radius of the cluster 
is imbedded. Peaks in the range of position angle ~ 100° to ~ 200° represent the tidal 
bridge-like feature appearing in the surface density maps. An extratidal extension feature 
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toward the Galactic center is correlated with peaks in the range of ~ 250° to ~ 50° on the 
azimuthal number density profiles. 



4.4. NGC 5466 



To examine the spatial configuration of stars around NGC 5466, we first constructed 
a star count map, as shown in the top-left panel of Figure [161 with a spatial bin size of 
2.88' X 2.88' for ~ 3° x 3° total field of view. We then smoothed the star count map with 
Gaussian kernel widths of 0.168°, 0.264°, and 0.312°, yielding smoothed surface density 
maps shown in Figure [161 Isodensity contour levels in Figure [16] correspond to values of 
the background level and la, 2cr, 3a, Aa, Qa, 10a, and 40cr above the background. Contour 
levels in the star count map correspond to those of a smoothed map with a Gaussian 
kernel width of 0.312°. The circle centered on the cluster indicates its tidal radius of 



Tt = 20.98' ( Lehmann &: Scholz 



19971). 



As can be seen in Figure [161 ^in overdensity feature above the background level 
apparently extends beyond the tidal radius of NGC 5466. Moreover, the extending 
extratidal overdensity feature seems to form two tails in the southeast and in the northwest 
directions. The tails stretch out symmetrically to both sides of the cluste r, tracing over 



2° o: 



fl2006ah and 



arc in the observed field. 



Grillmair &: Johnson 



lese are the same features reported by 



Belokurov et al 



(|2006[ ) in which the orientation of the tails is consistent 
with the cluster's orbit, as judged from the proper motion data. Note that the long arrow 
in Figure [161 indicates the direction of the projected proper m otion of NGC 5466 , i.e., 



fiaCos6 = —3.90 ± 1.00 mas yr ^, fis = 1.00 ± 1.00 mas yr ^ ( lOdenkirchen et al 



19971), 



while the short arrow represents the direction of the Galactic center. The tails also weakly 
show an S-shape feature of tidal stripping for cluster stars, as appeared for Palomar 5 



( lOdenkirchen et al. 



20031 ). There are clumpy overdensity features in the southern area of 
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the ob served field, which were also found in the surface density map of iGrillmair fc Johnson 



( 120061 ) ■ However, neither the cluster's orbit nor the direction to the Galactic center is likely 



to have a correlation with the clumpy structure. The appearance of such clumps can be 
caused either by the background galaxies or by the local substructures of the stars in the 
Galactic halo. 

The upper panel of Figure [17] shows the radial surface density profile of NGC 5466, 
for which we measured mean stellar number densities in concentric annuli. The solid curve 
characterizes the empirical representation of the data fitted by the King model of c = 1.2. 
While the measured surface densities have a nearly constant completeness value of ~ 98%, 
the very central region of the cluster inner than log{r') ~ 0.1 seems to show crowding 
effects. Although the inner part of the profile is well represented by the King model, the 
observed profile departs from the model prediction at the radius of ~ O.Sr^, showing the 
extended overdensity feature to be ~ 2rj. The radial surface density profile in this region 
is characterized by a power law with a slope 7 = —2.44 ± 0.13. The me asured slope is 



signifi cantly steeper than the cases of a constant mass- loss rate 7 = — 1 fjJohnston et al. 



I999I ) and of the other target clusters in this paper. This is certainly because N GC 5466 



shows a well-developed structure of extratidal tails. Indeed, 



Rockosi et al. 



(12002h found the 



annular-averaged surface densities of stars along the well-defined extrat idal tails of Pal 5 are 



also fi tted to a power law in radius with a significantly steeper slope. 



Lehmann fc Scholz 



(119971 ) also found an increased surface density in the outer region of the radial density 
profile for NCG 5466, which may be an indication of a tidal tail caused by unbound stars 
around the cluster. The departure of the radial profile from the prediction of the King 
model is more obviously presented in the residual plot of Figure [T71 The lower panel of 
Figure [17] shows the radial surface density profiles of eight angular sections in the observed 
field, as mentioned in Figure IU The break at ~ 0.5r( and the increased density features 
out to ~ 2rt are shown in profiles for all eight angular sections, and are in good agreement 
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with the extratidal overdensity feature in the surface density maps of Figure [161 Also, the 
estimated mean surface densities (/i) of the radial range are relatively higher in angular 
sections 3, 4, and sections 7, 8 than in the other sections. This is clearly a consequence of 
the two tails in the southeast and in the northwest directions around the cluster, as shown 
in the Figure [T6l 

The orientation of the detected two tails is further analyzed in Figure [181 which shows 
annular (left panels) and radially cumulative (right panels) azimuthal number density 
profiles around NGC 5466. The area within the radial range of 15' ^ r < 75', i.e., 
0.7r+ < r < 3.6r+ of NGC 5466, was divided into five concentric annuli with radial bins of 
lO' (inner three) and 15' (outer two) widths. We then clockwise measured the azimuthal 
number density with a 10° angular bin. It is apparent from the inner two annuli of the 
Figure [TS] that there are two main overdensity ranges at the position angles ~ 120° to 
~ 210° and ~ 320° to ~ 30°, with peaks near 175° and 345°, respectively. The peaks seem 
to be slightly shifted to position angles ~ 140° and ~ 330° of the outer three annuli. These 
are well matched with the symmetric S-shape feature of the extratidal tails detected in the 
surface density maps of stars around NGC 5466, as shown in Figure [161 

5. Discussions and Summary 

Using homogeneous wide-field observations obtained with the Megacam mosaic imager 
on the 3.6 m CFHT, we have undertaken several analyses to identify spatial configurations 
of stars surrounding five metal-poor globular clusters in the Galactic halo, i.e., M15, M30, 
M53, NGC 5053, and NGC 5466. Because of the low metallicities of the clusters we expected 
to easily detect the extratidal features from background stars with the secured wide-field 
deep photometric data covering an approximate 3° x 3° field-of-view. Analyses of surface 
density maps revealed that all the clusters present features of extratidal overdensities 
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over their tidal radii in the form of extended tidal tails and halos around the clusters. 
The observed extratidal extensions likely trace the spatial orbits of the clusters as well as 
the effects of dynamical interaction with the Galaxy. Radial surface density profiles and 
azimuthal number density profiles confirmed the overdensity feature and the extratidal 
structure of tails and halos around the five metal-poor globular clusters. 

M15, which is placed at 10.4 kpc from the Galactic center, is o ne of the densest 



globu lar clusters in the Milky Way. It is a post core-co lla psed cluster (de Marchi &: Paresce 



19961 ) with a high concentration c = 2.5 (.Harris 



19961). 



Lauchner et al 



fl2006h pointed out 



that M15 is unlikely to have extratidal structures due to its high concentration. However, 
we found apparent extratidal features from the analyses of the spatial stellar distribution 
around M15. Indeed, surface density maps revealed features of tails in the directions of the 
Galactic center and anticenter and of tails in the direction of the cluster's orbit inferred 
from the proper motion. The radial surface density profiles of stars surrounding M15 have 
a break at ~ O.Tr^ and an exces s of surface den s ity ex tendi ng to ~ LSr^, which indic ate 



tidal tail features around M15. iGrillmair et al. 



(119951 ) and iLehmann fc Scholzl (119971 ) also 



found an indication for tidal tails in the radial stellar density profiles of M15. Radial surface 
density profiles in different directions and azimuthal number density profiles confirmed 
directions of the extended features of the tails around M15, which are seen in the surface 
density maps. From the results, we can conclude that the apparent tails around M15 are 
real tidal features of the cluster. No Abell clusters were found in the observed field, and 
the observed field was not polluted by a strong dust extinction (E{B — V) = 0.10). We 
note, however, that clumps in the observed field around M15 could be due to unidentified 
background galaxy clusters, although significa nt contamination by the c lustered background 



objects seems a priori unlikely in the field (cf. 



Odenkirchen et al 



200ll ) 



M30 is located on 7.1 kpc from the Galactic center. The cluster is one of the most 
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metal-poor globular cluste rs in the Mil ky Way with [Fe/H] = —2.12, and it has a very high 



concentration of c = 2.5 (IHarris 



19961 ). Extended extratidal features on the configuration 
of stars surrounding M30 have not yet been reported. Analyses of spatial surface density 
maps of stars in the vicinity of M30 revealed flocculent overdensity features with tails 
around the cluster that extend to about 2rt and beyond. The radial surface density profiles 
of stars have a break at ~ O.Gr^ for M30, and the surface density excess seems to extend 
to ~ 2rt. The apparent east- west extension of the extratidal tails is likely associated with 
the direction of the Galactic center and anticenter. This is consistent with the extratidal 
overdensity property appearing in the radial surface density profiles and the azimuthal 
number density profiles for stars surrounding M30. Although tails associated with the 
direction of proper motion do not clearly show in the spatial surface density maps, the 
radial density profiles in O.Gr^ < r < 2rt for M30 shows relatively shallower slopes with a 
power law in the direction of the cluster's spatial orbit. We note, however, the visibility of 
such tails around M30 depends on the orientation of the cluster's orbit to our line of sight, 
i.e., tails could not be observed if we were situated in the plane of a cluster's orbit. The 
background features of galaxies and dusts are not likely to disturb the observed flocculent 
overdensity features of tails around M30. Indeed, the observed field is polluted only by two 
Abell clusters. The dust extinction is low toward this cluster (E(i? ~ V) = 0.03) and the 
fluctuation of the dust emission is low as traced by the IRAS lOO/xm map. 

The clusters M53 and NGC 5053 are close to each other, with a projected angular 
distance of about 1°, so they constitute an unusual astronomical pair in the outer halo 
of the Milky Way. These two clusters are the most metal-poor halo globular clusters 



with similar metallicities, i.e., [Fe/H] = -1.99 for M53 and [Fe/H] = - 2.29 



5053 . Also they have the same relative age as M92 ( iHeasley &: Christian 



1991 



br NGC 



Rey et al. 



19981 ) and a nearly the same distance modulus. The diffuse cluster NGC 5053 has a tidal 



radius of about 13.67 with a very low concentration parameter c = 0.84, while M53 has 
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a rather higher concentration c = 1.78 and a tidal radius of 21.75' ( iHarria Il996l ). In 
this study, we identified extratidal features on the spatial configuration of selected stars 
in the vicinity of two clusters M53 and NGC 5053. The observed tails likely trace the 
clusters' orbits and the direction to the Galactic center. As new findings, we detected 
features of clumps and ripples surrounding the clusters, and a tidal bridge-like feature 
and an e nvelope structur e arou nd the multiple globular cluster system. Here, we also 



note that 



Lauchner et al 



(120061 ) reported a discovery of a 6° tidal stream extending from 
NGC 5053 in the southwest direction, using the SDSS data from which they removed stars 
within the tidal radius of M53 from the analysis. Unfortunately, it is difficult to confirm 
the extended stream in the surface density maps of Figure [11], because our observations do 
not contain photometric data in the southwest direction of the M53 and NGC 5053 system. 
Such spatial stellar configuration features might be a result of tidal interaction or merging 
between the two clusters. Indeed, the structures of clumps and ripples are p redicted from 



recen t numerical simulations for the interaction of two globular clusters (e.g., 



Miocchi et al. 



20061 ). Furthermore, features such as a tidal bridge and an envelope in a multiple globular 
cluster system are also predict ed from the numerical simulations o f the possible merging of 
a set of globul ar clusters (e.g. , 



point of view. 



Capuzzo Dolcetta fc Miocchi 



20081 ). From an observational 



van den Berghl (1l996l ) suggested that a merging of clusters would occur more 



preferentially in dwarf g alaxies with low velocity dispersions of the cluster system in the 



range of 6.5~11 kms ^ (iMateo 



19961 ) rather than in the Milky Way. 



To test the contamination of background galaxies in the field of M53 and NGC 5053, 
we analyzed the spatial distribution of non-pointlike SDSS sources. For the extracted SDSS 
sources, we transformed the SDSS g'r'i' magnitude system into Ci and C2 indices, using 
the Equation (1). Sources falling into the masks in (02,^ ) and (ci,z') planes, which were 
used to select cluster member stars, were applied to the sample of background galaxies. In 
Figure [191 we show a smoothed surface density map of the selected background galaxies. 
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overlaid with the surface density contours of cluster stars of Figure [TTl Filled triangles 
indicate two Abell clusters in the observed field. As can be seen in Figure [121 we found no 
strong correlation between the pattern of the surface density variation of the background 
galaxies and the observed spatial configuration of stars around M53 and NGC 5053. Also 
note that the dust extinction is low toward the clusters (i.e., E{B — V) = 0.02 for M53 and 
E{B — V) = 0.04 for NGC 5053), and the variation in dust emission through the observed 
field is low as traced by the IRAS 100 /xm map. 

According to photometric results in this paper, it could be suggestive that the clusters 
M53 and NGC 5053 may be interacting or have interacted in the past. In the kinematic 
point of view, however, it would not be airtight that they originated with a same common 
dwarf spheroidal parent. The radial velocity of M53 is —79.1 ±4.1 km s~^, but it is 



+44.0 ± 0.4 km s ^ for NGC 5053 (jHarrid 119961 ). In this case, the clusters might have a 



common space orbit to the similar direction, only if the tangential velocities are in same 
directions with exorbitantly larger values than their absolute values of radial velocities. 
The tangential velocity of M53 is estimated to be /i = 43.1 ± 84.0 km s^^ by using values 
of the proper motion, i.e., /IqCos^ = 0.50 ± 1.00 mas yr~^ and fig = —0.10 ± 1.00 mas 



yr ^ (jPinescu et al. 



19991 ). Even if we are taking into account of the large measurement 
error, the calculated tangential velocity of M53 is merely comparable to the absolute value 
of the radial velocity. Unfortunately, we can not estimate the tangential velocity of NGC 
5053 because the proper motion has not yet been reported. Thus, the current kinematic 
information about radial velocities and inaccurate tangential motions of M53 and NGC 
5053 makes it still hard to arrange if the clusters came from a common dwarf and were on 



similar paths through the Gala ctic ha 



M53 and 16.19 for NGC 5053 ([Harris 



o. Als o, their distance moduli, (m — M)^ = 16.31 for 



19961 ). are still separated by ~ 1.4 kpc, a significant 



distance in the halo of the Milky Way co mpared with their tidal radii 112.7 pc for M53 



and 66.8 pc for NGC 5053 ([Harris 



1996[). Nevertheless, it should be noted that radial 
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velocities of M53 and NGC 5053 are commensurate wit 
in the Magellanic plane (see Figure 4c in 



Yoon fc Lee 



1 tho se of the other clusters aligned 



2002I ). Indeed, radial velocities of 



the globular clusters in the plane show a sinusoidal pattern of Keplerian motion as they are 
plotted against the orbital longitude measured along the Magellanic plane. This is possibly 
indicative of a common origin of the clusters. 

The cluster NGC 5466 is at a 16.2 kpc distance from the Galactic center. It is one 



of the lo w-met a 
c = 1.2 (IHarris 



icity 



19961) 



clu sters with [Fe/H] = —2.22 , and it has a low concentration with 



Lehmann fc Scholzl (119971 ) found evide nce of tidal tails around 



NGC 5466 in a radial density profile of the cluster. Recently, 



Belokurov et al. 



fl2006bh 



observed tidal tails of NGC 5466, stretching ab out 2° in apposite directions. Using SDSS 



photometric data. 



Grillmair fc Johnson 



( I2OO6I ) found evidence for a much larger extension 



of the tidal tails, with the leading arm over ~ 30° and t he trailing arm extending at least 



15°. From numerical simulations. 



Fellhauer et al 



(I2OO7I ) reproduced the tidal tails around 



cluster NGC 5466. In our study, analyses for spatial configuration of stars around NGC 
5466 have indicated an overdensity feature extending outside the tidal radius of the cluster. 
In the spatial surface density maps in Figure [161 we have found an S-shape extratidal 
feature of tails stretching out symmetrically over ~ 2° in the observed field. The observed 
radial density profiles of NGC 5466, as shown in Figure [TTl show a departure from the 
King model, with a break at ~ 0.5rt extending the overdensity feature out to ~ 2rt of 
the cluster. Azimuthal number density profiles in Figure [18] confirmed the orientation of 
the two tails around NGC 5466, which is likely consistent with the cluster's orbit inferred 
from the cluster's proper motion. The results are reconfirm ation of the tida l tails o f NGC 



5466 stretching as much as 45° on the sky as reported by 



Belokurov et al. 



Grillmair fc JohnsonI (l2006l ). 



(12006bh and 



Similar to the case of M53 and NGC 5053, we checked contamination of background 
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galaxies in the field of NGC 5466 by analyzing the spatial distribution of non-pointlike 
SDSS sources. Figure [2U] shows a smoothed surface density map of the selected background 
galaxies overlaid with the surface density contours of cluster stars from Figure Uni Filled 
triangles indicate eight Abell clusters in the observed field. As shown in Figure [201 the 
pattern of density variations in the selected galaxy sample does not seem to be correlated 
with the location of the tidal tails around NGC 5466. However, the spatial configuration 
of the selected background galaxies in the southern area of the observed field is likely 
correlated with the clumpy overdensity feature of the surface density maps of selected 
cluster member stars of Figure [161 The observed field is not polluted by strong dust 
extinction (E(5 - V) = 0.00). 

The other two low-metallicity clusters M68 and M92 in the Magellanic plane, for which 
we did not secure the wide-field deep photometric data from the present observations, 
are located on 10.1 kpc and 9.6 kpc from the Galactic center. Indications of extratidal 
extension s of stars ar o und t he clusters incidently have been overlooked in a few previous 
studies. 



Testa et al. 



(l20o"oi used the DPOSS plates to show an circular struct ure of 



extratidal halo around M92 extending out to ~ 30' from the cluster center. Also, 



Lee et al. 



(I2OO3I ) used CFHT mosaic CCD imaging to discover a presence of extratidal materials in 



the radial dens it y pro file of M92 and a weak tidal halo around the cluster. In addition. 



Grillmair et al. 



(I1995I ) found an indication for extratidal overdensity feature around M68 
using photographic plates taken at the UKST. Nevertheless, there have been little evidences 
for extratidal tails toward any preferred directions on the spatial stellar distributions 
around both of the clusters because of the lower spatial resolutions and the brighter limiting 
magnitudes of the photographic and CCD photometry. Thus we note the need for deeper 
wide-field CCD photometry to investigate spatial structure of the extended extratidal tails 
around both of the clusters. 
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The spatial configuration of stars tidally stripped from a parent stellar system of a 
present globular cluster contains direct information of the dynamical evolution that the 
cluster has experienced. Along with internal dynamical processes, external perturbations, 



such as disk an d bulge shocks 



evolution (e.g.. 



jy the Galaxy, could accelerate the cluster's dynamical 



Spitzer fc Hart 



1971 



Gnedin fc Ostriker 



19971 ). With respect to external 



dynamical evolution, three clusters, M15, NGC 5053 and NGC 5466 among our five targets 
might have experienced strong interactions with the Galaxy in the form of disk and bulge 
shocks, as indicated by the r atios of destruction rates with Vtot/vevap = 3.3, 130, and 160, 



respectively, for each cluster ( IGnedin fc Ostriker 



19971 ). Here, Vtot and Vevap are the total 



destruction rate and the evaporation rate per Hubble time, respectively. Furthermore, 



LeonetaL 



(120001 ) suggested that tidal tails usually stretch into the Galactic center, as the 
effect of both bulge and disk shock or the bulge shock to the cluster is much stronger than 
disk shock for the cluster. This implies that tidal tails stretching into the Galactic center 
and anticenter around the three clusters M15, NGC 5053 and NGC 5466 might be a result 
of both the bulge and disk shock of the Galaxy. Two other clusters, M30 and M53, are not 



expected to suffer strong gravitational shocks by the Galaxy, with Vtot/v, 



evap 



1.0. This 



would be a reason why the tidal tail feature of M30 is not well defined in the direction of 
the Galactic center. Also, the tidal bridge-like feature and an envelope structure around 
M53 may not be the results of gravitational shocks of the Galaxy but the interaction with 
the companion cluster NGC 5053. In addition, apparent extratidal extensions around the 
other two globular clusters M68 and M92 with Vtot/vevap = 1-37 and 1.00, respectively, may 
not be the results of suffering strong gravitational shocks by the Galaxy. In spite of such a 
dependency on the Galactic potential to the dynamical structure of the observed tidal tails, 
we found here that the feature of extratidal extensions around the target clusters is unlikely 
to have any distinctive correlation with the Galactocentric distance of each cluster. 



On the other hands, many numerical simulations of tidal tails performed over the 
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past decade (i.e., 



Montouri et al 



Combes et al. 



1999 



Yim fc Lee 



2002 



Capuzzo Dolcetta et al 



2005 



20071 ) have suggested that the stars dissolved from a globular cluster by the 
Galactic tidal force slowly drift along the cluster's orbit due to the orbital motion and the 
coriolis acceleration, and then eventually form a dense tidal tail parallel to the cluster's 
orbit. In this paper, we found that all of the observed low-metallicity clusters have shown 
possible extratidal extensions to the direction of the clusters' orbits, which are inferred 
from their proper motions. Indeed, the five target globular clusters are assumed to lie on 
one spatial orbital pla ne in the outer Ga lactic halo, and their kinematic motions place also 



almost on this plane (lYoon fc Lee 



20021 ). which are considered as a re sult of accret i on of 



dwarf galaxies. Based on a comparison of population and kinematics. 



Yoon fc Led (120021 ) 



suggested that 



recent studies ( jPalma et al. 



jMC is the possib 



2002 



e parent galaxy of these halo globular clusters. How ever, 



Bellazzini et al. 



2003 



Martinez-Delgado et al. 



20M) have 



stated that M53, NGC 5053 and NGC 5466 are possible former members of the Sgr galaxy. 
To examine a possible association of the clusters with the Sgr orbit or Magellanic debris 
plane, we mark the three clusters on the planes formed by the rectangular Galactocentric 
coordin ates with X, Y, an d Z in kiloparsecs, as shown in Figure Uli The Magellanic debris 



plane (lYoon fc Lee 



19981 ) are also overlaid in 



20021 ) and the Sgr orbit (jibata fc Lewis! 
Figure EH We also plotted the space motion vectors of two clusters M53 and NGC 5466 
with the uncertainties in directions, in order to see their velocity coherence to the two 
orbital planes. Space motion vector of NGC 5053 is not plotted, since the proper motion 
of the cluster is not available. As shown in Figure [21], M53 is likely to move in a direction 
along the Magellanic plane. The space motion of M53 is perpendicular to the Sgr orbit in 
X — Y plane, and it is directionally opposite in X — Z plane to the orbital direction of the 
Sgr dwarf. The space motion vector of NGC 5466 in X — Z plane seems to be parallel to the 
direction of the Sgr orbit, but it is not conclusive that NGC 5466 is moving along the Sgr 



orbit. The space motion of NGC 5466 in X — Y plane is against the Sgr orbit, and it seems 
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to m ove in a direction a long to the Magellanic plane within the uncertainty order of ~ 2.5a 
(e.g., 



Yoon&Lee 



2002). In addition, NGC 5466 is closer in space to the Magellanic plane 
than to the Sgr orbit. NGC 5053 lies remarkably closer to the Magellanic plane than to the 
Sgr orbit, while it has no information about the space motion. These lead us to conclude 
that the three clusters M53, NGC 5053, and NGC 5466, are potentially associated with the 
Magellanic plane rather than the Sgr orbit. 

Under these observational conditions, tidal tails for the five target globular clusters 
are expected to have a common feature of kinematic orbital properties, which would 
provide information about their parent dwarf galaxies. Although the spatial extensions 
and morphological shapes exhibited different regimes of developments for the observed 
extratidal tails, all the observed target clusters presented evidences of tidal tails tracing 
their orbits, as inferred from the projected directions of their proper motions. Thus, it 
appears that the tidal tails to the cluster's orbit for five target metal poor halo clusters 
are indeed evidences for an accreted origin from satellite system. Also, these observational 
results of extratidal tail features around the five target clusters provide further evidences of 



the merging scenario of the formation for the Ga. 



constraint on the potential in the Galactic halo (jOdenkirchen et al. 



actic halo, and in additi on an important 



200% . 



Data for more accurate proper motions of target clusters and radial velocities of stars 
in the observed extratidal structures are essential to verify their association with a stream 
in the Galactic halo. Furthermore, larger area coverage with accurate deep photometric 
data will also enable us to better constrain the orbits of the target clusters and cluster's 
member stars in the extratidal structures. With information about the spatial positions, 
radial velocities and proper motions of stars in the tidal tails, we could integrate their orbits 
backwards in the Galactic potential and eventually ascertain the origin of the globular 
clusters in the Galactic halo. 
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Fig. 1.— The Digitized Sky Survey image for the field of M53 and NGC 5053. Boxes 
represent the observed nine Megacam fields. Circles indicate the tidal radii for each cluster. 
North is up and east to the left. 
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Fig. 2. — The (c2, i') and (ci, i') color- magnitude diagrams of stars in the fields of five target 
globular clusters. From left to right, the first panel is the {c2,i') plane for stars in the 
empirically determined central regions of each cluster, the second to fourth panels are the 
(ci,i') planes for stars in the clusters central area, in the assigned backgournd fields, and in 
the total observed field, respectively. The fines in the (c2,i') plane indicate 2ac2 limits at i' 
magnitudes due to photometric error. The lines in the (ci, i') planes are the selected subgrid 
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Fig. 3. — Spatial distribution of all stars measured in the vicinity of five globular clusters. 
Circles centered on each cluster indicate the cluster's central areas in which the number ratio 
of stars enclosed in 2a di^') and those outside of 2ac^{i') in the (ci,i') plane in Figure [2] has 
a maximum value. The other four circles in each field indicate the background areas for 
guaging the field star contamination. 
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Fig. 4. — The schematic reseau plot used for the star counts. The radial surface densities 
were measured in concentric annuli. Each annunus was divided into eight angular sections 
(S1~S8) to derive surface density profiles for a different direction. Azimuthal number den- 
sities were determined in an assigned annulus with respect to a position angle measured 
clockwise from the east principal axis with a 10° unit. 
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Fig. 5. — From top-left to bottom-right, the raw star count map around M15 and the surface 
density maps smoothed by Gaussian kernel values of 0.090°, 0.135°, and 0.195°, overlaid with 
isodensity contour levels of 2a, 2.5a, 3a, Aa, 8a, and 30cr above the background level. The 
circle indicates a tidal radius of M15. The short and long arrows represent the direction of 
the Galactic center and of the cluster's proper motion, respectively. 
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Fig. 6. — Upper : Radial surface density profile of M15 with a theoretical King model of 
c = 2.5 (solid curve). Residuals, log(0/C), after subtracting the values predicted by the King 
model from the observed surface densities are also plotted. Radial completeness, the ratio of 
the number of recovered artificial stars and the total number of stars, and its trend (dotted 
line) are added. The profile within the detected overdensity region, ~ 0.7rj < r <~ 1.5rf, 
is represented by a power-law, i.e., a dotted straight line in logarithmic scale with a slope of 
7 = —1.59 ± 0.21. The mean number density of stars in the overdensity region is estimated 
as /i = 0.943 per arcminute square. The dashed horizontal line indicates the background 
density level, and an arrow represents the tidal radius of M15. Lower : Radial surface density 
profiles of eight different angular sections (S1~S8). For clarity, the surface densities from SI 
to S8 are given zero-point offsets with values of const. = 8.4, 7.2, 6.0, 4.8, 3.6, 2.4, 1.2, 0.0. 
Dotted lines indicate the power-law fit to the overdensity region of ~ O.Tr^ < r <~ 1.5rt 
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Fig. 7. — The annular (left panels) and radially cumulative (right panels) azimuthal number 
density profiles of M15 with an angular bin size of 10°. The radial bins are set to be lO' for 
the inner three annuli and 15' for the outer two annuli. The horizontal dashed lines are the 
determined backgound density level and dashed dotted hues are the level of lu deviation 
above background level. 
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Fig. 8. — The raw star count map around M30 and the surface density contour maps 
smoothed by Gaussian kernel values of 0.084°, 0.147°, and 0.210°. The contour levels are 
the background level, 0.5cr, 1.5o", 2.5cr, 3.5cr, and 8a above the background level. The circle 
indicates a tidal radius of M30. The short and long arrows represent the direction of the 
Galactic center and of the cluster's proper motion, respectively. 
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Fig. 9. — Upper : Radial surface density profile of M30 with a theoretical King model of 
c = 2.5 (solid curve), residuals after subtracting the values predicted by the King model, and 
radial completeness and its trend (dotted line). The profile within the detected overdensity 
region, ~ 0.7rt < r <~ 2rt, is represented by a power-law, i.e., a dotted straight line in 
logarithmic scale with a slope of 7 = —1.41 ± 0.21. The mean number density of stars in the 
overdensity region is estimated as yU = 0.410 per arcminute square. The dashed horizontal 
line indicates the background density level, and an arrow represents the tidal radius of M30. 
Lower : Radial surface density profiles to eight different angular sections (S1~S8). For 
clarity, the surface densities from SI to S8 are given zero-point offsets with values of const. = 
10.5, 9.0, 7.5, 6.1, 4.7, 3.3, 1.5, 0.0. Dotted lines indicate the power-law fit to the overdensity 
region with 7 as a slope, fi is the mean density in the overdensity region. 
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Fig. 10. — The annular (left panels) and radially cumulative (right panels) azimuthal number 
density profiles of M30 with an angular bin size of 10°. The radial bins are set to be lO' for 
the inner three annuli and 15' for the outer two annuli. The horizontal dashed lines are the 
determined backgound density level and dashed dotted hues are the level of la deviation 
above background level. 
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Fig. 11. — The raw star count map around M53 and NGC 5053 and the surface density 
contour maps smoothed by Gaussian kernel values of 0.125°, 0.175°, and 0.240°. The contour 
levels are the backgound level, 2cr, 4o", 8a, and 20cr above the background level. The circles 
indicate tidal radii of M53 and NGC 5053. The long arrow from M53 indicates the direction 
of the cluster's proper motion. The short arrows represent the direction of the Galactic 
center. 
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Fig. 12. — Upper : Radial surface density profile of M53 with a theoretical King model of 
c = 1.78 (solid curve), residuals after subtracting the values predicted by the King model, and 
radial completeness and its trend (dotted line). The profile within the detected overdensity 
region, ~ CTr^ < r <~ l.Gr^, is represented by a power-law, i.e., a dotted straight line in 
logarithmic scale with a slope of 7 = —1.58 ± 0.19. The mean number density of stars in the 
overdensity region is estimated as /i = 0.302 per arcminute square. The dashed horizontal 
line indicates the background density level, and an arrow represents the tidal radius of M53. 
Lower : Radial surface density profiles of six different angular sections. For clarity, the 
surface densities from SI, S2, S3, S4, S7 and S8 are given zero-point offsets with values of 
const. = 8.1, 6.7, 5.2, 3.9, 1.3, 0.0. Dotted lines indicate the power-law fit to the overdensity 
region with 7 as a slope, fj, is the mean density in the overdensity region. 
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Fig. 13. — The annular (left panels) and radially cumulative (right panels) azimuthal number 
density profiles of M53 with an angular bin size of 10°. The radial bins are set to be lO' for 
the inner three annuli and 15' for the outer two annuli. The blank bins correspond to the 
areas without photometric data and the area within the tidal radius of the neighbor cluster 
NGC 5053. The horizontal dashed lines are the determined backgound density level and 
dashed dotted lines are the level of Icr deviation above background level. 
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Fig. 14. — Upper : Radial surface density profile of NGC 5053 with a theoretical King 
model of c = 0.84 (solid curve), residuals after subtracting the values predicted by the King 
model, and radial completeness and its trend (dotted line). The profile within the detected 
overdensity region, ~ O.Tr^ < r <~ 2.5rt, is represented by a power-law, i.e., a dotted 
straight line in logarithmic scale with a slope of 7 = —0.62 ± 0.15. The mean number 
density of stars in the overdensity region is estimated as /i = 0.235 per arcminute square. 
The dashed horizontal line indicates the background density level, and an arrow represents 
the tidal radius of NGC 5053. Lower : Radial surface density profiles of eight different 
angular sections (S1~S8). For clarity, the surface densities from SI to SB are given zero- 
point offsets with values of const. = 12.6, 10.8, 9.0, 7.2, 5.4, 3.6, 1.8, 0.0. Dotted lines indicate 
the power-law fit to the overdensity region with 7 as a slope, fj, is the mean density in the 
overdensity region. 
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Fig. 15. — The annular (left panels) and radially cumulative (right panels) azimuthal number 
density profiles of NGC 5053 with an angular bin size of 10°. The radial bins are set to be 
10 for the inner four annuli and 15 for the outer two annuli. The blank bins correspond 
to the areas without photometric data and the area within the tidal radius of the neighbor 
cluster M53. The horizontal dashed lines are the determined backgound density level and 
dashed dotted lines are the level of Icr deviation above background level. 
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Fig. 16. — The raw star count map around NGC 5466 and the surface density contour maps 
smoothed by Gaussian kernel values of 0.168°, 0.264°, and 0.312°. The contour levels are 
the backgound level, la, 2cr, 3a, 4cr, 6cr, lOa, and 40(7 above the background level. The circle 



indicates a tidal radius of NGC 5466 (ILehmann &: Schob 



19971 ). The short and long arrows 



represent the direction of the Galactic center and the cluster's proper motion, respectively. 
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Fig. 17. — Upper : Radial surface density profile of NGC 5466 with a theoretical King 
model of c = 1.20 (solid curve), residuals after subtracting the values predicted by the King 
model, and radial completeness and its trend (dotted line). The profile within the detected 
overdensity region, ~ O.Sr^ < r <~ 2rt, is represented by a power-law, i.e., a dotted straight 
line in logarithmic scale with a slope of 7 = —2.44 ± 0.13. The mean number density 
of stars in the overdensity region is estimated as /i = 0.375 per arcminute square. The 
dashed horizontal line indicates the background density level, and an arrow represents the 
tidal radius of NGC 5466. Lower : Radial surface density profiles of eight different angular 
sections (S1~S8). For clarity, the surface densities from SI to S8 are given zero-point offsets 
with values of const. = 14.0, 12.0, 10.0, 8.0, 6.0, 4.0, 2.0, 0.0. Dotted lines indicate the power- 
law fit to the overdensity region with 7 as a slope, fj, is the mean density in the overdensity 
region. 
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Fig. 18. — The annular (left panels) and radially cumulative (right panels) azimuthal number 
density profiles of NGC 5466 with an angular bin size of 10°. The radial bins are set to be 
lO' for the inner three annuli and 15' for the outer two annuli. The horizontal dashed lines 
are the determined backgound density level and dashed dotted lines are the la deviation 
above background level. 
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Fig. 19. — The gray scale map for the surface number density of background galaxies around 
M53 and NGC 5053. Background galaxies were selected by using non-pointlike SDSS sources. 
Filled triangles indicate Abell clusters in the observed field. The overlaid contour levels show 
the spatial configuration of selected stars surrounding two clusters, M53 and NGC 5053. 
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Fig. 20. — The gray scale map for the surface number density of background galaxies around 
NGC 5466. Background galaxies were selected by using non-pointlike SDSS sources. Filled 
triangles indicate eight Abell clusters in the observed field. The overlaid contour levels show 
the spatial configuration of selected stars surrounding NGC 5466. 
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Fig. 21. — The spatial distribution of seven low- metalli city globular clusters in the Galacto- 
centric coordinate. The Galactic plane is drawn in X — Z plane, and the Sun (indicated by 
the circle with a dot) is at (-8,0,0) kpc. The filled circles are five target clusters in this study, 
while two open circles represent M68 and M92. Space motion vectors of M53 and NGC 5466 
are indicated by thick arrows with the uncertainties. The filled s quare marks the S agittarius 
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Table 1. Basic data for five target globular clusters 



Target 


a 


5 


1 


b 


Rsun 


Rgc 


rc 


rt 


[Fc/H] 




(J2000) 


(J2000) 


(deg) 


(deg) 


(kpc) 


(kpc) 


(') 


(') 




M15 


21:29:58.3 


12:10:01 


65.01 


-27.31 


10.3 


10.4 


0.07 


21.5 


-2.26 


M30 


21:40:22.0 


-23:10:45 


27.18 


-46.83 


8.0 


7.1 


0.06 


18.34 


-2.12 


M53 


13:12:55.3 


18:10:09 


332.96 


79.76 


17.8 


18.3 


0.36 


21.75 


-1.99 


NGC 5053 


13:16:27.0 


17:41:53 


335.69 


78.94 


16.4 


16.9 


1.98 


13.67 


-2.29 


NGC 5466 


14:05:27.3 


28:32:04 


42.15 


73.59 


15.9 


16.2 


1.20 


20.98 


-2.22 



and Rgc represent distances from the Sun and the G alactic center, respectively, rc 

1 1996 ) except for core 



Note. — Rsun 

and rt are the core radius and tidal radius of a cluster. Data are from 



Harrii 



radius and tidal radius of NGC 5466, which are from 



Lehmann fc Scholj 119971) 



Table 2. Journal of Observations 



Target 


Filter 


Exp. Time 


FWHM 


Number of 






(s) 


{") 


fields 


M15 


9 


135 


0.8 


9 




r 


170 


0.8 


9 




i 


260 


0.8 


9 


M30 


9 


110 


0.7 


8 




r 


125 


0.7 


8 




i 


160 


0.6 


8 


M53 / NGC 5053 


9 


90 


1.0 


9 




r 


180 


1.0 


9 




i 


480 


0.9 


9 


NGC 5466 


9 


90 


0.9 


9 




r 


180 


0.8 


9 




i 


480 


0.7 


9 
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Table 3. The coefficients a and h for a rotational transformation of the {g —r',r' — i) 

plane to {ci,C2) plane, for each cluster. 



Target a b 



M15 0.905 0.425 

M30 0.909 0.417 

M53 / NGC 5053 0.909 0.416 

NGC 5466 0.912 0.411 



